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Abstract. The well-known analytical solution of Burgers’ equation is extended to 
curvilinear coordinate systems in three dimensions by a method that is much simpler 
and more suitable to practical applications than that previously used [22], The results 
obtained are applied to incompressible flow with cylindrical symmetry, and also to 
the decay of an initially linearly increasing wind. 


1. Introduction. Burgers’ equation is a well-known example of a nonlinear partial 
differential equation whose solution can be constructed from a linear partial differen- 
tial equation. It is, to the best of our knowledge, the only such example. Much of the 
interest in this equation arises because it is a very simple form of the Navier-Stokes 
equation in the one-dimensional, Cartesian, time-dependent, compressible, viscous 
limit [17]. An early study derived two steady-state solutions [1], and then the equa- 
tion was used in studying the decay of free turbulence [4, 5], The importance of 
the equation is due to the nonlinear term, udu/dx , which allows the calculation of 
the modification of the velocity due to the exchange of momentum between a vari- 
ety of different length scales. Burgers’ equation is the simplest type that shows the 
complicated interplay between the nonlinear steepening and the diffusion of a wave. 

The first full solution of the one-dimensional Burgers’ equation was found inde- 
pendently by both Cole [7] and Hopf [12], Recently, the solution has been extended 
to n dimensions by using group action on coset bundles [22]. The basic idea of this 
approach is to determine coset spaces for a chosen group and for a chain of closed 
subgroups, and then to construct a bundle out of these spaces. Next, the action of 
the group on the bundle is prescribed and a class of tensor-valued functions called 
invariants on the coset bundle is introduced. Because these group invariants take the 
form of differential equations, the method has been used to obtain a tensor version 
of the n-dimensional Burgers’ equation and to show that its solutions can be con- 
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structed from a linear tensor diffusion equation. It has also been demonstrated that 
the procedure gives a nonlinear tensor equation for a generalized Burgers’ equation 
with a higher than second order derivative. In this paper, we examine extensions 
o the basic class of solutions of the vector Burgers’ equation, which we discovered 
before we became aware of the work of Wolf et al„ to more than one dimension and 
to curvilinear coordinate systems. The main advantage of our approach is that it is 
extremely simple from a mathematical point of view and that the results we obtain 
have the potential to be used in a number of applications (see Section III for two 
simple applications). It must be added, however, that the solutions obtained in this 
paper represent a subclass of more general solutions to the tensor Burgers’ equation 
discussed by Wolf et al. [22], In what follows, we recapitulate part of the historical 
derivation presented first by Cole [7] and Hopf [12]. This derivation is central to our 
results in Sections II and III. 

The equation to be solved is 

du du d 2 u 

dI + u d^ = l 'J7-' (1) 

where v is a physical constant (viscosity). A successful transformation uses 

d(j> 

U ~fo' <t> = 4>{x , t). (2) 

Substituting into Eq. ( 1 ) and integrating term-by-term with respect to x , 

dcf> 1 (d<j>\ 2 d 2 <t> 

IS?) < 3 > 

where the function of integration is omitted. Because Eq. (3) is invariant under the 
transformation 



( 8 ) 
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Quite an interesting literature has developed around applications of the solutions 
in pure mathematics [20], magnetohydrodynamics [9], astrophysics [14, 16], and 
cosmology [10, 15]. Lighthill [18] extended the basic transformation to the case 
where the x-coordinate is measured in the frame of reference moving in the same 
direction as the wave at the undisturbed speed of sound. A valuable book on Burgers 
equation and related topics was recently published but is, unfortunately, already out 
of print [11]. 

In the following, we formally present the derivation of the solutions to the vector 
Burgers’ equation by using a generalized Cole-Hopf transformation. 

2. The vector extension of Burgers’ equation. We begin with the vector equivalent 
of Burgers’ equation 

^ + w • Vm = vV 2 u. (9) 

at 

In order to evaluate the vector Laplacian we must use 

V 2 m = V(V • u) - V x V x m, (10a) 

and for the inertial term 

2 

u ■ Vu — ^2 — m x V x ii. (10b) 

Then Eq. (9) greatly simplifies as long as 

u x V x u = L'V x V x u. (10c) 

This appears to be an equation that represents quite general flows, but the only 
solutions are irrotational. This can be seen by taking the divergence of both sides of 
Eq. (10c) and using 

V • (i*x B) = B • V x A- A- V x B. (10d) 

In what follows we will assume that the flow field is irrotational since this is necessary 
to simplify Eq. (9) and to derive the extension of the Cole-Hopf transformation. The 
assumption is satisfied trivially for the one-dimensional Cartesian problem 

V x m = 0, (Ha) 


so that 

w = V<£. (Hb) 

Eq. ( 1 lb) is the extension of Eq. (2) to three dimensions and, together with Eqs. (10), 
allows the simplification of Eq. (9) to 


V 


'd± , 
dt 2 



= 0 . 


(12) 


The meaning of Eq. (12) is that the quantity enclosed by brackets must be a function 
of time, only, because its gradient is zero in all space so that 


d± , (V0) 2 

dt 2 


-vS/ 2 4> = E{t). 


(13) 
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The form of the unknown function E(t) is irrelevant for our discussion since it 
cannot affect the velocity. We simply define a new potential 


4>i =<t>- j E{t)dt, 


where the gradients of <f> and <j> x are identical. 

Before introducing the Cole-Hopf assumption, it is useful to state the following 
vector identity: 


V« = f|v 0 , 


(14a) 


which is true for any orthogonal, curvilinear coordinate system. Taking the diver- 
gence of both sides with a = </> and using Eq. (14a) again with a — jjj yields 


v >= ?-^(V0 ) 2 + ||v 2 0. 


sr 


de 


(14b) 


We will also need the diffusion equation 


dd 2 


(15) 


Now we apply Eqs. (5), (14b), and (15) to Eq. (13) to generate Eq. (7a), again so 
that 


h = - t V0. (, 6) 

This clearly shows that the Cole-Hopf transformation goes through for quite general 
coordinate systems. Cole [7] originally discussed the vector solution for the case of 
Cartesian coordinates, and he indicated that it applied to some special solutions in 
higher dimensions. 

The usefulness of Eq. ( 16) can be shown by exploiting the simple symmetry of the 
following example. 


3. Cylindrical symmetry. Using only a radial velocity. Burgers’ equation in cylin- 
drical coordinates with axisymmetry may be written as 


du du 
~di + U Tr =V 


1 d / du\ u 

[rdr \~dr / ' (17 > 

We initially noticed that (17) could be written as 

du du d(\d \ 

J7 + u a? ! “'a-r{7a? (r ">)- < 18 > 

Equation (18) and the related integrals of (18) are what led us to believe that the 
Cole-Hopf transformation was quite general in nature. The particular form of ( 1 6) 
is 

-2udd 

U = ~d~~dr (> 9 ) 

where the diffusion equation is 


dQ_ _ v_d_ / dd\ 
dt r dr \ dr ) ’ 


(20) 
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A simple example of the use of this method is to use the initial-value solution of 
the diffusion equation (a Dirac delta at / = r = 0 ) 

6 = J-e~^. ( 21 ) 

2 vt 

Then the solution of (17) is found from (19), 

u = -. (22) 

t 

This is the cylindrical version of Burgers’ Cartesian solution [4], u = x/t , which is the 
inviscid solution of (1) (solution 2.1, figure 5 in Benton [3]). Because (22) satisfies 
(17) in the limit as v -*■ 0 (which is also the case of zero acceleration for fluid 
particles), this solution is the inviscid limit of the solution of Burgers’ equation in 
cylindrical coordinates. Some solutions of (17), such as example A below, approach 
the sawtooth limit as / -► oo ; but this is not surprising since it has been noted [2, 
1 3] that the sawtooth shape of dissipation layers is of fundamental importance in 
Burgers’ representation of turbulence. 

It is extremely useful to derive solutions of the diffusion equation that are based 
on physical boundary conditions on the velocity field. To this end, Eq. (19) may be 
formally integrated to 


6(r , t) = k(t) exp J u(w , /) < 


where 

k(t) = 0(0,t). (24) 

The velocity is to be initially specified over a given range so that the initial value of 
6 is given by 


6(r, 0) = 0 o (r) = k 0 exp j u Q { 


w) dui 


for r 0 < r < R . 

The general solution of the diffusion equation, (20), with azimuthal symmetry, is 

1191 , , 

»(',')=[ I H^r )J < /kr')J,fkr)c krdkdr. (26) 

Jo Jo 

The integral over k can be done once and for all since it is a special case of Weber’s 
second exponential integral as discussed in Watson [21]: 

j" e\p(-i'tk 2 )J 0 {kr)J 0 (kr l )k dk = ^ exp h fe) ( 27 > 


where the modified Bessel function of order zero may be written as 

I 0 {ar) = J 0 {iar), 

which is a real function. Equation (26) may now be written as 


6{r,t) 


■ exp - 


l l±\r'dr'. 
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4. Two simple examples. We choose to work two relatively simple examples so 
as not to obscure the simple nature of the transform. It is surprisingly easy to pose 
apparently straight-forward problems that cannot be readily carried to completion 
because of extremely difficult quadratures. We will treat the more difficult case of 
the decay of an initially constant velocity wind in a future paper. 

The reader may have noticed by now that Burgers’ equation does not specify the 
density distribution except in a few limits besides the trivial one of constant density 
(as in the next example). 

A. Incompressible flow. The first case is that of incompressible flow where we take 
a source of water emanating from the origin of coordinates and spreading out in a 
circularly symmetric pattern. The initial flow is assumed to exist from very close to 
r = 0 out to R , which is assumed to be larger than any other length scales in the 
problem. We then take the limit as R -* oo , in the same spirit as other studies of 
this initial-value problem on infinite intervals [2]. The speed must decrease as r~' 
to conserve mass in an initially steady flow. In sum, 


u(r, 0) = k 


r n <r<R, 


so that (25) integrates to 


where 




Clearly, u 0 may be large at the small value r 0 in such a way that the product u 0 r {] 
is well defined. Equation (29) can now be integrated using HankeFs generalization 
of (27) as discussed in Watson [21], so that 

*-VST(i-f)c (33) 

where (§) n is Pochhammer notation, 

(i).-!G + 0(I*’)-(5—') <34. 

and 

(f) 0 = 1 <35) 

The full nonlinear solution can now be derived from (19) and expressed in terms 
of the confluent hypergeometric function as 


or more clearly as 


M(r , t ) = a ' ^(a/2+E2,-r 2 /(4w)) 

2t M(a/2, 1 , -r 2 /(4w)) 


E °° { ~ \ n ( r \ n 

n=l\2>nT^?\ 4l/ ,l 

M(r, t) = -122 

t OO 1 (Q \ ( r \n 

*~~*n=0 (,j| v 2 / 4 


(36b) 
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Isochrones of Nonlinear Solution «(r, /) 



Fig. 1. The nonlinear evolution of a wind with an initial spatial 
dependence of \jr . The isochrones show the decay of the wind, 
and. especially, the diffusion at the leading edge where the wind was 
turned off as this discontinuity propagated outward. 

We see that the solution can be written as the product of the inviscid solution times 
an infinite sum of powers of the similarity variable r /(4 vt) . We note the nonlinear 
effect on higher powers of the similarity variable through the occurrence of the factor 
n in the numerator. The transient nonlinearity steepens but diffusion dominates; 
note that the transient is not by itself a solution of Burgers’ equation since we do not 
have the superposition theorem of linear theory. 

Figure 1 was created by nondimensionalizing Eq. (17) in the usual fashion, i.e., 
r —* r/L, t —* vt/l 2 , u -* u/U . Choosing U = m 0 and L = r Q , then a is 0.5. This 
ratio of series is initially divergent and required 100 terms in double precision on a 
Sun SPARCstation. The solutions are presented as isochrones, showing the evolution 
of the nonlinear solution towards a steady-state. The boundary where the wind was 
turned off moves outward while diffusion smooths out the discontinuity. The wind 
initially drops off as 1 jr so that the peak in the profile approaches the starting radius 
for the source as t — * 0 . All curves approach the steady-state envelope of 1 /r while 
diffusion flattens the leading edge. Asymptotically the velocity is zero everywhere. 

B. The decay of an initially linearly increasing wind. Example two is a wind em- 
anating from the origin, again, but now the velocity linearly increases in magnitude 
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out to some radius, r Q , chosen to be larger than any other length scale in the prob- 
lem. We had in mind the initial theoretical treatments of radiatively-driven winds [6] 
from early-type stars, but our description is too simple to be applied to this problem. 
This is because of our assumption of cylindrical rather than spherical symmetry, and 
also because other effects are now known to be important in these winds [8], 

The driving force is imagined to turn off at t = 0 , and then we calculate the decay 
of the flow. In sum 

u{r, 0) = 0<r<r ; t = 0, (37a) 

'o l 0 

where 


is half the time for the wind to flow from the source to the initial outer boundary. 
Both r 0 and u 0 are large, but t 0 is not necessarily big. Then 


Proceeding as before, 


V 


-hr 


b = 


4i >r n 


(38) 

(39) 


6 = 


l+'/' o eXP 

u{r, t) = 


1 


4 vt 1 + tjt 

r 


t + u 


We note that the solution approaches the Burgers’ sawtooth shape as t 


(40) 

(41) 


OC . 


5. Conclusions. We have given a relatively simple proof that Burgers’ equation 
can be solved for curvilinear coordinate systems in three dimensions by assuming 
irrotational flow and using Eq. ( 1 6) with d derived from solutions of the diffusion 
equation. We solved two examples in cylindrical coordinates with circular symmetry 
to show how to implement the generalized solution. Our goal has been to extend the 
tools of workers who model the complicated interplay between nonlinear steepening 
and viscous diffusion in waves. 
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